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Abstract 

The possibility of detection of some features of high energy particle interactions with 
detectors placed at medium depths underground through studies on high energy muons is 
studied. These muons carry information about the early interactions occurring during the 
development of the hadron cascade near the top of the atmosphere. They might reveal the 
effects resulting from creation of quark gluon plasma in interactions of ultra high energy 
cosmic ray iron nuclei with air nuclei. 
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1. Introduction 



The cosmic ray composition near the knee region (5 x 10 14 eV up to almost 10 17 eV) 
is a challenging problem [[IJ. Recent measurements by KASCADE have suggested 
an increased proportion of heavier nuclei in cosmic ray flux at these energies. However, 
these results are still far from conclusive as the results derived from electron, muon and 
hadron data differ and the composition measurements from Cerenkov arrays show even 
decreasing trend ||, Q. Although more precise determination of the chemical composition 
is not yet feasible, the presence of Fe 56 nuclei at ~ 10 — 20% level in the primary flux 
at ~ 10 15 eV seems quite plausible. This result presents an interesting opportunity to 
observe iron nitrogen and iron oxygen nucleus-nucleus interactions at energies beyond 
the possibilities of terrestrial accelerators. Though both N 14 and O 16 are light nuclei^, 
the number of participating nucleons in head-on collisions with iron nuclei may go upto 
several tens. Thus, given the very high energy per nucleon available with the incoming iron 
nuclei, the formation of quark-gluon plasma seems to be a likely possibility in a fraction 
of the primary interactions. Although the baryon density is only half of the density in 
Pb-Pb interactions (due to the interaction with lighter nuclei), the energy density e can 
reach sufficiently high value. This will be shown in section ^] The conditions in which 
such interactions are taking place are unique. However they provide an opportunity to 
investigate some characteristics of quark-gluon plasma (QGP) which are not available 
presently in terrestrial experiments. The interactions of iron nuclei take place at altitudes 
higher than 20 km and secondary charged pions have large probability to decay into muons 
due to the thin air at these altitudes. While the decay products of neutral pions undergo 
further interactions leading to the development and absorption of the electron-photon 
casacde in the middle and lower atmosphere, the muons, in some sense, decouple from the 
shower development and are subjected only to energy loss through ionization and multiple 
scattering. Moreover, the large Lorentz factor helps in keeping the secondary particles 
collimated within a narrow cone. As a result, there is a possibility to observe differences 
in the characteristics of the muon component in some extensive air showers from what 
one would expect on the basis of standard models of high energy interactions. 

The fact that the hard component of muons in extensive air showers carries informa- 
tion on the pion production in the initial stages of the shower development is of great 
importance. Hence, aside from QGP, it is also interesting to study the sensitivity of ener- 
getic muons to the details of the primary interaction dynamics. In this regard the QGP 
effects observed in the muon component of air showers will also indicate the extent of the 
above mentioned sensitivity. Of course it is necessary to determine in more quantitative 
way the notion of 'hard component' in terms of energy and the methodology to observe it. 
One simple method is to place a large area high resolution muon detectors underground at 
such depth that the overburden will filter out the less energetic muons and preferentially 
provide the information on the initial interaction. 

Due to the absence of a comprehensive theory on the formation and evolution of QGP 
and also due to the lack of solid experimental evidence there are many phenomenological 
models with quite different and sometimes even opposite predictions. The same is also 
true for results on secondary particles multiplicity. There are predictions of a decrease 



' In principle interactions with Ar i0 may also happen but these are suppressed by a factor larger than 
~ 100 according to CORSIKA simulations. 
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of multiplicity 0] in high energy nuclear collisions. On the other hand, many statisti- 
cal thermal models predict increase of multiplicity J7|. Anyway the secondary particle 
multiplicity is an important variable and a change in the interaction dynamics may be 
observable through change in the characteristics of particle production, depending on the 
size of the effect and the sensitivity of measurements. The aim of this paper is to study 
such phenomena through detailed simulations and to learn about conditions necessary 
to observe at least some of the features which would reveal the formation of QGP in 
cosmic ray interactions in the atmosphere. To achieve this, a very simple model of QGP 
formation is presented in section ^] This model is then used to simulate the formation 
of QGP and the results are compared with the outcome of programs currently used for 
simulations of extensive air showers. The analysis of simulated data and comparisons are 
presented in section |3] The conclusions are discussed in the final section. 



2. Quark Gluon Plasma Simulation 

The standard package used for simulation of cosmic ray interactions in the atmosphere is 
the program CORSIKA |§ supplemented by several models of high energy interactions 
of elementary particles. The properties of these models have been extensively studied 
in ||, [jj~0| . Good agreement with experimental data of KASCADE has been achieved using 



QGSJET ||rj and VENUS [LJ models with QGSJET giving somewhat better results. 
VENUS predicts slightly larger number of muons per shower than QGSJET at energies 
in the range, 10 14 eV - 10 16 eV. As this feature mi ght blur or hide the muon component 
contributed by the QGP, the VENUS model has been used for comparisons while searching 
for the observable effects due to the QGP. Another reason for this choice is the availability 
of good documentation for VENUS and the fact that QGP simulations supplemented by 
VENUS simulations were possible without any intervention with the VENUS code. 

The key factor in QGP modelling is the number of nucleons, N int , participating in 
the first interaction of an iron nucleus with the atmospheric nucleus. If Ni nt is larger 
than a certain limit, N^ot, nucleons are supposed to melt down into QGP. The values for 
the fraction of melted nucleons used in simulations were 10 %, 20 % and 40 % and the 
corresponding threshold values of N int were chosen to be 13, 15 and 19 nucleons. We shall 
denote this fraction as TZh = N hot /N int . Thus, in the case of TZh = 0.1, only two nucleons 
can melt. For values of 0.2 and 0.4, this minimum is 3 and 8 nucleons respectively. The 
aim was to keep the minimum Ni nt leading to QGP production to be above 10 nucleons. 
In the process of simulation, Nhot nucleons are taken aside and rest of the particles par- 
ticipates further in collision simulated with the VENUS code. The total energy of N hot 
nucleons, denoted E hot , is used to produce again N hot nucleons, but this time with the 
kinetic energy generated according to the Boltzmann distribution and isotropically in the 
overall centre of mass system. The remaining Ehot is used for production of ir^ and 7r°, 
again with Boltzmann distributed kinetic energy and isotropically in the centre of mass 
system. The pion production is stopped when the E hot is reduced to 2.5 GeV . Afterwards 
this remnant cluster is passed to CORSIKA for further processing. For each of the three 
values of TZh considered, the pions were generated with Boltzmann distribution parameters 
corresponding to the mean pion momenta < p 7T > values of 60 MeV, 120 MeV, 300 MeV 
and 1 GeV. The wide range of the values of the fraction, TZh, together with a range of 
possible values of the mean momentum of generated pions are intended to account for 
uncertainties in the description of the QGP freeze out. The properties of eventual QGP 
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< P 7T > 
MeV 


TZh ratio 


0.1 


0.2 


0.4 


60 


1492 (880) 


2743 (1717) 


6764 (3802) 


120 


1389 (791) 


2610 (1634) 


6332 (3568) 


300 


1074 (680) 


2165 (1398) 


4918 (2614) 


1000 


674 (393) 


1260 (825) 


3093 (1808) 



Table 1: Mean numbers of pions (both and it ) produced in QGP. Figures given in 
the parentheses represent the R.M.S. values. 

formed in light nuclei interactions are also unknown. For example, the temperature might 
not be high enough to produce cc or even ss quark pairs. Also the isotropy of pions ra- 
diated during the freeze out is not assured and the mean momentum along the shower 
axis may differ from momenta in the transverse plane. However, it turns out that the 
results are not sensitive to this. The results of simulations with < >= 300 MeV and 
< > = 30 MeV do not differ significantly from those with overall < p 71 >= 300 MeV. 
Hence, for the sake of simplicity, the pion production is considered to be isotropic. The 
mean number of all kinds of pions produced under different conditions in this simple 
approach are given in the Table |1| The number of pions produced per nucleon-nucleon in- 
teraction and its dependence on the primary energy are shown in Fig. [I]. The two extreme 
cases - one with maximum pion production, i.e. 40 % of iVj ni nucleons melting down and 
pions evaporating with mean momentum of 60MeV (henceforth we shall abbreviate this 
as QGP(0A, 0.06)), and the other with minimum pion production (QGP(0.1, 1.0)) are 
shown. All the data sets under consideration here were generated from 2000 iron nucleus 
initiated showers with the primary energy within the range 10 15 eV — 2. 10 17 eV (the 
upper limit is the one recommended as a limit of reliability of the VENUS program ||). 
The spectral index of the energy distribution was —2.7 M, The primary zenith angle was 
fixed to be degrees while the azimuth angle and the altitude of the initial interaction 
were left free. The mean primary energy of 2.2 10 15 eV is relevant to all results presented 
in the paper if not stated otherwise. 

The chosen primary energy interval should not be taken as a prediction for QGP for- 
mation at any specific energy. The primary energy range was determined essentially by 
the range of validity of the VENUS model and by flux considerations, i.e., the experiment 
with realistic area should be able to observe the effects. This point will be further dis- 
cussed in section fl] The kinematical conditions and the dynamics of parton interactions 
as described by VENUS allow the QGP formation at these energies. One can estimate 
the energy density e using the Bjorken's expression 0] 

AE/Ay 
€ ~ hR 2 t ' 

where R is the radius of the smaller nucleus, r is the formation time (l-r2 fm/c) and E is 
the c.m.s. energy of final state particles in the rapidity unit interval near y crns = 0. With 
the transition temperature in a two-flavour QGP in the range, 150 MeV < T c < 200 MeV, 
the corresponding critical energy densities are 0.8 GeV fm~ 3 < e c < 2.5 GeVfm~ 3 . The 
VENUS simulations give for Fe — N interactions at Ei nc = 10 15 eV e = 3.3 GeV fm~ 3 
with impact parameter b < 2 fm and r = 1 fm/c (e = 3.5 GeVfm~ 3 for b < 1 fm). 
At energy E inc = 10 16 eV, the energy densities reach the values e = 4.9 GeV fm~ 3 for 
b < 2 fm and e = 5.4 GeV fm~ 3 for b < 1 fm. Thus even with some uncertainty in the 
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Figure 1: Number of pions (both 7r ± and 7r°) produced per NN interaction. 
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value of r, the creation of QGP in the iron-air interactions at the "knee" energies is not 
excluded, at least according to VENUS. The lower part of the primary energy interval is 
already accessible at RHIC. However, the collider experiments explore mainly the central 
rapidity region while the cosmic ray experiments are sensitive to particle production in 
the very forward region. Thus the two approaches are in a sense complementary. 



3. Quark Gluon Plasma Effects 

The main detectable quantities used to describe extensive air showers are the electro- 
magnetic and muon sizes and corresponding lateral distribution functions (LDF). These 
quantities are related as follows, e.g. for the muon size iV M , 

= f 2 2Trp^r)dr, 
J n 

where p M (r) is the density of particles in the detection plane and r is the distance from 
the shower axis. 

About one third of pions originating from QGP are tt°. Their decay products un- 
dergo further interactions in the atmosphere and they will disperse among other particles 
of electromagnetic nature. As the electromagnetic size is typically much larger than the 
muon size, the QGP effects should be easier to detect in the muon component of extensive 
air showers. However, the muon sizes of showers at the ground level do not reveal any 
differences. The truncated muon sizes Nj[ are compared in Fig. ^. The truncation corre- 
sponds more to the reality of the experimental conditions where the measurements near 
the shower core are hampered by systematics caused by punch-through and depend on the 
details of detector design. The kinetic energy cutoff used here for pions is E kin > 300MeV. 
The results obtained from simulations with QGSJET as well as VENUS but without QGP 
formation are compared in the figure with simulations of QGP supplemented by VENUS 
in the case of the most abundant pion production, QGP(0A, 0.06). 

Even the more detailed lateral distribution functions shown in Fig. ^| are almost iden- 
tical and do not reveal any features which would allow us to observe effects corresponding 
to QGP formation. The events with QGP have slightly larger LDF values in the region 
50 4- 100 meters than the VENUS simulation. However, this difference is not measurable 
and the differences between QGSJET and VENUS are larger than any QGP effects. The 
differences for events with other values of IZh and < p 7T > are similar or even smaller. 

As a consequence of kinematics, the muons originating from QGP formed in the initial 
interaction of iron nucleus with air have mostly higher energies than the muons produced 
lower down in the atmosphere. Thus detection of hard muons in cores of extensive air 
showers should reveal differences between cases with and without QGP production. To 
estimate this effect, the data are further analysed with detectors at a medium depth 
underground. Our somewhat simplified virtual experiment is placed at a depth of 110 m 
with 100 m thick overburden made of standard rock (A = 22, Z = 11 and p = 2.65g/cm 3 ). 
The underground cavern with the 10 m high ceiling is supposed to have detection area 
of 100 x 100 m 2 . The tracking of particles generated by CORSIKA through the rock 
was done using GEANT |13|]. The rock filters out the muons with momenta smaller than 



53 GeV. 

The muons penetrating underground have an average momentum around 280 GeV 
but the distribution has a very long tail. To illustrate the contribution of QGP originated 
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Figure 2: Muon size truncated between 30 and 230 meters from the shower centre at the 
ground level. 
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Figure 3: Lateral distribution of muon densities on the ground. 
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model 


model only 

[GeV] ' 


QGP + VENUS 
[GeV] 


QGSJET 

VENUS 

QGP(0.4,0.06) 

QGP(0.4,0.12) 

QGP(0.4,0.3) 

QGP(0.4,1.0) 


287.5 (265.5) 

287.2 (259.8) 
99.9 (31.1) 

105.0 (33.6) 

130.3 (41.4) 
234.5 (73.7) 


277.8 (225.8) 

277.9 (226.8) 
279.3 (232.6) 
282.2 (244.8) 



Table 2: Mean momentum of muons detected underground in iron initiated showers. The 
R.M.S. values are given in the parentheses. 

muons relative to standard muon production, Fig. |] shows the momentum distribution 
of muons produced in a typical event with E inc = 2.25 10 15 eV, which is close to the 
average primary energy, and with the height of the primary interaction at 33.1 kilometres, 
corresponding to the maximum of the distribution. In this particular event Ni nt = 36, 
which is higher than the average. The reason for this selection was to obtain at IZh = 0.4 
high QGP muon multiplicity. The superimposed contributions from QGP(0A, 0.06) and 
QGP(0A, 1.0) show strong dependence of QGP muon momenta on the initial conditions 
of pion production. The underground detection area is sensitive only to a small fraction of 
muons coming from QGP pion decays as can be seen from Fig. |5|. The overall momentum 
distribution is changed only slightly as may be seen from the Table 0. The values given in 
this table are the averages for the interval 0-=-30 TeV to cut away the extreme fluctuations. 
Close values for QGP(0A, 0.06) and QGP(0A, 0.12) are given by the momentum cutoff - in 
the case of QGP(0A, 0.06) one can observe underground only the tail of the distribution. 

The underground muon size (necessarily truncated by the design of experimental 
area) reveals notable differences already in the case of unmodified models (Table |3|). The 
value used to define the tail of the underground muon size distribution, 990 muons, has 
been chosen so that in the case of pure VENUS simulations 400 events pass the limit 
and the statistical error is only 5 %. The effects of QGP are clearly detectable when the 
fraction of melted nucleons becomes larger. The muon size at medium depth underground 
also shows differences between VENUS and QGSJET in the case of proton showers. The 
muon LDF obtained from VENUS have systematically higher values, both for iron nuclei 
and protons, than those generated by QGSJET (Fig. |). In all cases compared here, 
the statistical errors are smaller than the marker size. Apparently such an underground 
experiment can essentially help to tune the models which is an indispensable condition 
for any further studies on chemical composition of cosmic rays and search for QGP effects 
in cosmic ray showers. 

The underground muon LDF with QGP simulation differs from the reference LDF 
obtained from data simulated by VENUS code only. Fig. |7] a)-d) show the underground 
muon lateral distribution functions corresponding to different values of the ratio TZh di- 
vided by the reference LDF. Each picture summarises the distributions with the same 
value of < p n >. The functions correspond to the events with large underground muon 
size - iV" > 990, where the differences were most notable. The LDF of VENUS and 
QGSJET with the same condition, > 990, are mutually much less different than the 
LDF without any conditions in Fig. - within the statistical errors they practically co- 
incide over the whole range 4- 55 m. Hence the effects caused by QGP production are 
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Figure 4: Momentum distributions of muons from QGP and from standard production. 
Altitude of interaction is 33.1 km, E inc = 2.25 10 l5 eV, N int = 36 and Nhot = 14. 
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Figure 5: Picture la,b): QGP(0.4,1.0); 2a,b): QGP(0.4,0.06). 1,2a) - gray: all possible muons; hatched: muons from decays within 
7 km from initial interaction; black: detected muons. 1,2b) gray: all QGP pions; black: pions whose daughter muons are detected. 







fraction 


of 


events 


model 


average 


with N$ > 990 
(in %) 


VENUS p 


464 


6.3 


± 


0.6 


QGSJET p 


378 


3.7 


± 


0.4 


VENUS Fe 


807 


20.1 


± 


1.1 


QGSJET Fe 


709 


13.6 


± 


0.9 


QGP(0.1,0.06) 


849 


21.8 


± 


1.2 


QGP(0.1,0.12) 


843 


21.7 


± 


1.2 


QGP(0.1,0.3) 


829 


20.4 
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1.1 


QGP(0.1,1.0) 


819 


20.4 


± 


1.1 


QGP(0.2,0.06) 


868 


23.4 


± 


1.2 


QGP(0.2,0.12) 


871 


22.9 
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1.2 


QGP(0.2,0.3) 


867 


22.5 
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1.2 


QGP(0.2,1.0) 


809 


19.3 
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1.1 


QGP(0.4,0.06) 


967 


27.0 
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QGP(0.4,0.12) 


955 


26.5 
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1.3 


QGP(0.4,0.3) 


904 


24.3 
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1.2 


QGP(0.4,1.0) 


835 


22.0 
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1.2 



Table 3: Average underground muon size of events and the fraction of events with 
All QGP simulations were supplemented by VENUS. 
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Figure 6: Lateral distributions for muon densities underground produced by VENUS and 
QGSJET. 
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in this sense model independent. Fig. [7] a)-d) show that the effects of more abundant 
pion production caused by increasing TZh can be offset by increasing < p n > to a certain 
extent. Near the centre of the shower, both scenarios give similar results. But at greater 
distances from the shower core the LDF with QGP contribution decreases slower than the 
reference one and starts to differ. This 'point of departure' moves away from the shower 
centre with increasing < p 7T > and becomes more prominent with the increase of TZh- 

As we are not dealing with real data it is not considered necessary to increase the 
statistics to pinpoint exactly the values of TZh and < p n > at which the distinction between 
models with and without QGP can be made. The comparison of the LDF's shows that 
with increasing < p n > one has also to increase the distance upto which the LDF can be 
measured in order to be able to see possible QGP effects. At small TZh one has to increase 
the observation time or the detection area to achieve the same result. Thus the resolving 
power of the experiment is related to the dimensions of the detector available underground. 
Assuming the detection area to be 100 x 100 m 2 and 10 % iron nuclei and 90 % protons in 
the primary flux, one can get in 2.4 10 7 sec (75 % of a year) 15 000 iron initiated events 
and 34 000 proton initiated events, with > 990. These numbers correspond to showers 
initiated by particles with E inc > 10 15 eV. In the case of E inc > 10 16 eV the overall flux 
is smaller and one would get only about 200 iron events and 400 proton events, assuming 
the same fraction for events above 990 muons cut as at E inc > 10 15 eV. Actually the 
percentage of events above this cut increases with energy and the yield of observed events 
should be higher. 



4. Discussion 

Results obtained from simulations discussed above cannot be compared at present with 
real data as the muon detectors are either located at ground level where the high en- 
ergy component is entirely hidden in the overall muon flux or they are located so deep 
underground that the muon momentum cutoff is of the order of TeV or larger and only 
few muons penetrate to the detector ( |jT3}| , p5] , fT6|| , [ |T7|j , fllBf , fl~9||). The few experi- 



ments with medium overburden have either crude spatial resolution like the underwater 
experiments |20[ and |21J or they have relatively small detection area like [^2J (6x6 m 2 ). 



The present study is based on the assumption that the high energy interactions of 
cosmic rays in the atmosphere may lead to the production of QGP in the initial stage 
of the shower development. The inclusion of a very simple model of QGP production 
and freeze-out into high energy interactions of iron nuclei with air in cosmic ray showers 
has shown that under certain kinematical conditions one could observe the QGP induced 
effects in the high energy muon component of extensive air showers. These kinematical 
conditions are mainly the availability of sufficient energy E^ot of melted nucleons and 
not too high mean momentum < p* > of the evaporated pions. The model has been 
intentionally chosen so simple so that it can cover sufficiently large interval of these two 
quantities with minimum parameters. Moreover the simulations have shown that the 
underground data on muons are rather insensitive to the isotropy of initial QGP pion 
distribution (Fig. |5|). It is interesting that none of the simulated QGP systems increases 
the muon density in the very centre of the shower and this number is essentially determined 
by the unmodified high energy interaction model (VENUS in this case). 

Even in the situation when there is no apparent difference in the underground muon 
size or underground LDF, like in the case of QGP(0.1, 1.0), the QGP production can 
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Figure 7: Ratio of underground lateral muon densities with QGP simulation to simulations 
done by VENUS only (events with > 990). 



12 



> 20 r 

CD 

°- 18 E- 

1 fi L VENUS only 

cn 

^14 — 

CD 

£ 12 - 

CO 
D 

I 8 

E 6 



10 



4 
2 




1000 2000 3000 4000 5000 

underground muon size N u „ 



> 20 r ^ 

CD 

°- 18 E- 

16 L QGP(0. 1,1.0) + VENUS 

i_ -i / _ 

CD 1 ^ 

S 12 - 

CO 

S 10 

I 8 j- ; ■ 

§ 6 E- 





1000 2000 3000 4000 5000 

underground muon size N U M 



Figure 8: Underground muon size versus primary energy of iron nuclei. 
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still be detected provided there is an independent measurement of the shower energy, 
e.g. the muon size Nff on the surface. This quantity is proportional to the energy of 
the incident nucleus and the proportionality holds even underground (Fig. H). Since the 
creation of QGP distorts the energy — muon size relation, a comparison of the muon 
sizes on the surface and underground may indicate deviation from the standard scenario. 
Events under the diagonal of proportionality E inc ~ are what one would expect - as 
the abundant pion production in QGP mimics higher energy of the projectile nucleus. The 
events above the diagonal are more intriguing. In this case clearly the muons produced 
from QGP are not able to offset the loss which is due to several factors. One of them 
is that in normal interactions the leading particle effect constrains the lateral energy 
dissipation. No such effect has been included into QGP simulation and this combined 
with sometimes very high altitude of the primary interaction leads to the spraying of the 
QGP originated muons on a very large area. Moreover the pions from QGP mostly do 
not start particle cascades and their energy is not converted into secondary particles any 
more. Thus one can take roughly the muon multiplicity deficit in a standard interaction 
as N deff = k std ■ E hot - q std while the gain in multiplicity from QGP is N QGP = (k QGP ■ 
Ehot — Qqgp) " fdet, where f det is the detection efficiency factor which is quite small (see 
Fig. H), kg-cp and k std are the slope parameters and qQ G p and q st d are the intercepts. 
The net result Nqgp — Ndeff is apparently a linear function with positive intercept and 
negative slope. At smaller Ehot while the intercept prevails one gets higher multiplicity 
than normal and when higher energy portion is lost from standard production one gets 
events with lower muon multiplicities underground. 

Another factor contributing to the multiplicity deficit is the exemption of Nhot nucleons 
from further interactions - there are only 56 — Nhot nucleons with the energy Ei nc • (56 — 
Nhot)/ 56- The baryons evaporated from QGP cannot provide enough muons to balance 
the loss of forward production initiated by Nhot nucleons in standard interaction. This 
part of the deficit is quite difficult to express quantitatively in terms of linear relations. 
Otherwise it would be much easier to measure chemical composition of cosmic rays ! In 
fact, events with lower muon multiplicity underground could be mistaken as products of 
interactions of lighter primary nuclei. Therefore, the resolution in energy measurement 
both at surface and underground will be crucial to make the distinction. 

The measurements on energetic muons in cosmic ray showers at medium depth un- 
derground can reveal information on the initial stages of interactions which are otherwise 
inaccessible through measurements on the surface. Provided the QGP is formed in in- 
teractions of lighter nuclei, the number of produced muons will grow linearly with c.m.s. 
energy per nucleon. The ratio IZh should, in naive expectation, also grow with the energy. 
Hence the observable effects are likely to become more apparent at higher energies (in the 
simple model used here, IZh was considered to be energy independent). Thus, in principle, 
one should be able to study some of the scenarios of QGP deflagration or detonation at 
the freeze out described in |23| . 
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